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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 2861 

ANALYTICAL INVESTIGATION OF ICING LIMIT FOR DIAMOND- SHARED 
AIRFOIL IN TRANSONIC AND SUPERSONIC FLCW 


By Edmund E. Callaghan and John S. Serafini 


SUMMARY 

Calculations have been made for the icing limit of a diamond air- 
foil at zero angle of attack in terms of the stream Mach number, stream 
temperature, and pressure altitude. The icing limit is defined as a 
wetted-surface temperature of 32° F and is related to the stream condi- 
tions by the method of Hardy. 

The results show that the point most likely to ice on the airfoil 
lies immediately behind the shoulder and is subject, to possible icing 
at Mach numbers as high as 1.4. 


INTRODUCTION 

As the operational speed of aircraft is increased through the 
transonic region, the frictional heating available to prevent the forma- 
tion of ice on the aircraft becomes an important quantity. A study of 
the probable icing limits of a high-speed airfoil in transonic and super- 
sonic flow was made in terms of the pertinent variables of flight Mach 
number and free-stream conditions of pressure and temperature. The pro- 
cedure used to determine the icing limit is based on the method presented 
in reference 1 for calculating the surface temperature of an insulated 
body (a body within which heat is not conducted from one section to 
another) running fully wet in an air stream. The icing limit at any 
point on the airfoil was assumed to depend upon maintaining a tempera- 
ture of 32° F at that point. The free-stream static temperature corre- 
sponding to this icing limit was calculated for each point on the airfoil 
for particular values of free-stream Mach number and pressure. 

The results presented herein were calculated for a sy mme trical 
diamond airfoil at zero angle of attack for a range of airfoil- 
thickness ratios from 0.02 to 0.10, pressure altitude from sea level to 
45,000 feet, and free-stream static temperatures to -40° F. 
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SYMBOLS 

The following symbols are used in this report: 

Cp specific heat of air at constant pressure, Btu/(lh)(°F) 
e vapor pressure, lb/sq ft 

k e coefficient of evaporation 

coefficient of heat transfer 
L latent heat of vaporization, Btu/lb 

M Mach number 

molecular weight of air 
m e molecular weight of water vapor 

p static pressure, lb/sq ft 

r recovery factor 

« 

T static temperature, °R 

t/c airfoil thickness ratio, ratio' of airfoil maximum thickness to 
chord 

Tq c minimum free -stream static temperature corresponding to ice -free 
9 condition on surface as defined in equation (la), °R 

AT-^- temperature increment resulting from frictional heating in the 
boundary layer, °R 

X ratio of specific heats of air, 1-400 
Subscripts: 

0 free -stream static conditions 

1 local conditions at the edge of boundary layer 


s 


surface 
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AMLYSIS 

The method of reference 1 offers a convenient means for deter mi ning 
the local surface temperature on a wetted' tody as a function of the 
free -stream conditions of temperature, pressure, vapor pressure, and 
velocity, provided the flow field about the body is known. Ease of cal- 
culation was achieved by a number of simplifying assumptions. The ’ 
restrictions caused by thes°e assumptions will be discussed in relation 
to the general heat balance for an unheated body flying through icing 
conditions. 

Several investigators have made quite complete analyses of 'the heat 
balance for a body flying through a cloud (references 2 and 3) • The 
generalized heat balance for such conditions is well known and has been 
stated quite concisely by Messinger (reference 2) . Some simplification 
of the general equation may be obtained by considering the case of an 
insulated body subjected to sufficient frictional heating to maintain 
the surface ice free. In this case, the heat balance is as follows: 

(l) Heat due to the frictional, or viscous, effects 

plus 

(2) Heat due to the kinetic energy of the water droplets 

plus 

(3) Heat of fusion 
equals 

(4) Heat lost by convection 
plus 

(5) Heat for evaporation of water 
plus 

(6) Heat required to raise temperature of water droplets from stream 
temperature to surface temperature 

The result of Hardy's relation (reference l) is obtained by equat- 
ing the frictional term (l) to the sum of the convective term (4) and 
the evaporation term (5). For Hardy's relation to be valid, terms (2) 
and (3) must be small compared with term (l), and term (6) must be small 
compared with term (4) or (5); or the sum of terms (2) and (3) must be 
nearly equal to term (6). 
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The he&t of fusion term (3) offers considerable complication to the 
problem. It is conceivable that the surface film might consist of a 
mixture of water and ice particles without ice actually adhering to the 
surface. The heat of fusion would necessarily have to be withdrawn from 
a portion of the film for such particles to occur and, because of the 
lack of knowledge of the percentage of the film that might be frozen 
before visible icing occurs, it would appear advisable to define an ice- 
free surface as being fully wet 'at 32° F with no ice particles in the 
film. This definition for the icing limit eliminates term (3) from the 
general heat balance and Hardy's relation becomes applicable if term (2) 
is approximately equal to term (6) or both are negligible. Term (2) is 
directly proportional to the impingement rate and the square of the 
flight speed; term (6) is directly proportional to the impingement rate 
and the temperature difference between the stream and the surface. If 
the impingement rate is small, as for the case of low liquid-water 
contents, then both terms (2) and (6) are negligible. For calculations 
of the icing limit, both terms are of the same order of magnitude an d 
because of the nature of the problem, the temperature difference between 
the surface and the stream increases with increasing flight speed. Both 
terms therefore increase with increasing Mach number; large values for 
either term are only encountered at large values of liquid-water content 
and flight speed.. \ Such a combination is highly improbable in actual 
flight since high speeds usually occur at high altitudes, and high alti- 
tudes imply small liquids water contents. 

The relation of Hardy applies in particular to a surface barely 
wetted with a very thin film free of ice particles but is probably 
nearly correct for the whole-range of practical interest for high-speed 
flight. 


Method of Calculation 

Wetted-surface temperature . - The equation for the temperature of 
a wetted surface in an air stream developed in reference 1 was written 
as 


T 1 - T s + ^1 " 


k e m e / e s 

Vk \Pl - e s 



( 1 ) 


where the subscript s refers to conditions at the surface and sub- 
script 1 to conditions at the edge of the boundary layer. From the 
energy relation for a viscous fluid 

h + + 11 m i 2 ] 

Substituting this relation in equation (l) and solving, for T s yields 
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Tg = T-l 1 + Mp 2 


kgiiigL / e s 
^•h m a c p \Pl “ e s 



Because interest in this problem is fixe'd on a definition of flight 
circumstances that provide locarl surface temperatures of 32° F, the 
terms T s , L, and e s are constants. Hardy observed (reference 1) that 
for the range of temperatures near 32° F the ratio of the evaporation 
coefficient k e to the heat -transfer coefficient k^ was very nearly 
1. A solution of equation (l) can.be obtained for, any one of the four 
variables ( e-^, p^, M-j_, and T^) provided the other three are known. 

From purely physical reasoning, however, an additional restriction exists. 
This restriction results from the fact that the air in a cloud is fully 
saturated at the static or free -stream conditions. For each value of 
stream-static temperature, the free-stream vapor pressure is therefore 
assumed a constant and equal to the saturated vapor pressure. If it is 
assumed that the flow about the body, outside the boundary layer, is 
accomplished with no change in phase, that is, no condensation or evap- 
oration, then Dalton's law of partial pressures applies and 

fl = £l 

e o Po 


If a reasonable value of the recovery factor r is assumed, say 
0.88, then equation (l) may be rewritten in the form shown below 


T-i (l + 0.176 M-i ) 

492 = ; — T, 

1 0,c 


0,c 


- 2776.5 


12.75' 


_(| p o- 12 * 75 ) 


(Po - e o) 


(la) 


which is quite convenient for purposes of calculation because 
T l/ T 0,c ( 1 + °- 176 Mp 2 ) and Pp/p 0 can be readily calculated for each 
point on a diamond airfoil. 


Because eo is a function of To, c > the solution of equation (la) 
must be solved by trial arid error. In an effort to reduce the laborious - 
ness of such calculations, a. number of charts similar to that shown in 
figure l'were made. For each pressure altitude considered, the free- 
stream static temperature corresponding to an ice -free condition T n „ 

v/ y C 

was plotted as a function of the parameter Tp/To^c (l + 0.176 Mp 2 ) by 
means of equation (la) for constant values of the pressure ratio Pp/Pg* 
For each value of Pp/p 0 and T-l/Tq^ (l + 0.176 M-} 2 ) a free-stream 

static temperature which corresponds to a surface temperature of 32° F 
and therefore represents the minimum free-stream temperature for an ice- 
free surface T 0 ^ c is obtained. The lower limit of free-stream static 
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temperature was considered to "be -40° F for the purposes of these cal- 
culations since it has been shown (references 4_and 5) that supercooled 
water droplets are not likely to be found at temperatures below this 
value. 

Flow-field calculations . - The flow field about a diamond or double- 
wedge profile can be characterized in general by four regions; a subsonic 
region, a subsonic region with sonic velocity at the shoulder, a mixed 
sub sonic -supersonic region with a detached bow wave, and a supersonic 
region with attached shock waves. The regions of most interest are in 
the currently important ranges of transonic and supersonic flight an d 
the results presented herein will be limited thereto. 

A recent analysis is given in references 6 and 7 for the sub- 
sonic region with a sonic line at the shoulder. This solution is for 
a thin, symmetrical wedge at zero angle of attack followed by a straight 
section. The results in this region presented herein are therefore 
limited to the front half of the diamond airfoil. 

The mixed subsonic- supersonic region with detached bow wave has been 
solved in references 8 and 9. This solution applies to a thin, symmet- 
rical diamond airfoil at zero angle of attack for Mach numbers from 1 to 
the Mach number at which the shock attaches. The flow in this region is 
subsonic over the front half' of the airfoil increasing to sonic velocity 
at the shoulder. The flow then expands supersonically around the corner, 
and supersonic flow exists, over the back half of; the airfoil to the 
trailing edge. 

The flow about a diamond airfoil in purely supersonic flow has been 
known for many years. The results used herein are those of reference 10 
because of their ease in use and applicability. The flow in this region 
is supersonic over the front half of the airfoil with an expansion 
region at the shoulder and a consequent higher Mach number on the back 
surface of the airfoil, the velocities and pressures being constant over 
each surface for a particular value of airfoil thickness and stream Mach 
number. 

The calculation of the chordwise pressure and Mach number distri- 
butions over a symmetrical-diamond airfoil of specific thickness ratio 
t/c at zero angle of attack through the entire transonic- supersonic range 
can therefore be accomplished by combining the results of references 6 
to 10. 

The analyses of references 7 to 9 were made using a speed function’ f 
and a generalized pressure coefficient C p . The results of these analyses 
are presented graphically by plots of the chordwise distribution of | 
and C p for constant values of the free -stream speed function |q. The 
pressure ratio P,/Pq is related to the generalized pressure coefficient 

Cp by the equation 
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£l 

PO 


1 + 



2(r + i)- 3 


The Mach number, either local or free stream, is related to the 
speed function £ hy the equation- 


M = 

The temperature ratio Tj/T 0 is obtained from the energy relations 

for a compressible fluid and is given in terms of Mach number by the 
equation 



T x 1 + Mq 2 



An example of a typical calculation of Tq c at a particular 

chordwise position is given for the 25-percent chord station and a value 
of the free-stream speed function of 0.703- Values of £ of -0.740 

and Cp of 2.890 are obtained from figures 3 and 5, respectively, of 
* reference 8. If t/c of 0.06 is assumed, then Mq = 1-092, Mp = 0-8926, 
Pl/P0 = Tp/To = 1.068, and Hj/Tq (l + 0.176 M] 2 ) = 1.218. At a 

pressure altitude of 15,000 feet, a T 0 of 422.5° R is obtained from 
figure 1. 3 

For the supersonic case with an attached shock (reference 10), the 
values of P]/pq and M-^/Mq can be obtained directly from the charts of 

the reference report and, hence, the calculation of the parameters 
required for obtaining T Q ^ C from figure 1 is simple. 


RESULTS AND DISCUSSION 

For each value of airfoil thickness ratio and free-stream Mach 
number assumed, the values of the parameters T^/Tq (l + 0-176 Mp 2 ) and 

Pl/P 0 can be calculated over the whole range of transonic and supersonic 

speeds. At each pressure altitude, the free-stream static temperature 
T 0,c corresponding to a wetted surface temperature of 32° F (ice free) 

can then be obtained from charts similar to figure 1. 
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This value of T 0 ^ c corresponds to the lowest free-stream static 

temperature Tq which will give ice-free conditions at the particular 
point on the airfoil under consideration for the given values of free- 
stream Mach number and pressure altitude. 

A plot of Tq^ c as a function of the chordwise distance is shown 

in figure 2 for several free-stream Mach numbers for a 4 -percent-thick 
symmetrical-diamond airfoil at zero angle of attack and a pressure alti- 
tude of 30,000 feet. 

A study of figure 2 shows that the general trend of the curves ‘is 
very similar for Mach numbers from 0.7593 to 1.105, that is, from the 
subsonic region to the supersonic region with attached shock. For each 
curve the minimum T q^ c increases slowly from the leading edge to the 

midchord. The sudden jump in the free-stream static temperature required 
for icing protection at the midchord results from the expansion around 
the corner and the consequent decrease in local static pressure and 
increase in local velocity. The increase in velocity and the decrease 
in pressure enhance the evaporation, hence, the increase in the required 
To,c- Water impinging on the front half of the airfoil might therefore 
run back along the surface, spill around the shoulder, and freeze. 
Downstream of the shoulder, the free-stream temperature required for 
protection decreases slowly to the trailing edge. 

When the shock wave is attached and the flow is everywhere super- 
sonic (Mq = 1.2, or greater) ,. the free-stream temperature for protection 
is a constant for the whole front surface and the whole back surface. 

The free-stream temperature required for protection of the back half is, 
however, greater for the reasons previously cited. 

This and similar figures for other altitudes and airfoil thickness 
ratios show that the critical points for these airfoils are immediately 
ahead of the shoulder for the front half of the airfoil and immediately 
behind the shoulder for the back half of the airfoil. The icing limit 
for these two critical points is shown in figure 3(a) where the minimum 
free-stream temperature corresponding to an ice-free condition is shown 
as a function of the free-stream Mach number for pressure altitudes of 
sea level and 30,000 feet for a 6 -percent -thick airfoil. The value of 
T 0 , c decreases rapidly with increasing Mach number; that is, the tem- 
perature region of possible icing hazard becomes much smaller with 
increasing Mach number. The free-stream static temperatures required 
for protection apparently are higher at the higher pressure altitudes as 
would be expected since, with other conditions equal, the decreased 
pressure enhances the evaporation rate. At a pressure altitude of 
30,000 feet the curve shown for the critical point behind the shoulder 
indicates the possibility of icing at a Mach number approaching 1.4. 
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The effect of altitude is shown clearly in figure 3(b) for free-stream 
Mach numbers of 0.67, 1.0, and 1.3 for the critical point ahead of the 
shoulder on the front surface and for Mach numbers of 1.0 and 1.3 for 
the point immediately behind the shoulder. At a particular value of 
Mach number the curves for both critical points remain nearly parallel 
over the altitude range. For Mach numbers of 1.0 or less the free- 
stream temperature required for protection increases with increasing 
altitude with a definite decrease in the slope of the curves in the 
higher altitude range. For the higher Mach numbers the slope of the 
curves increase over the whole range from-sea level to 45,000 feet. 

The effect of airfoil thickness on the icing limit for the two 
critical points on the airfoil is shown in figure. 3(c) where the free- 
stream temperature required for protection is plotted as a function of 
the airfoil thickness ratio. ’ The icing limit at a pressure altitude of 
45,000 feet is shown for the point ahead of the shoulder for free-stream 
Mach numbers. of 0.875, 1.00, and 1.30 and for the point immediately 
behind the shoulder for Mach numbers of 1.0 and 1.3. For the point 
ahead of the shoulder, the temperature required for protection remains 
constant or decreases with increasing thickness depending on the Mach 
number. For the point immediately behind the shoulder, the temperature 
required for protection increases with increasing thickness. 

An increase in Tq^ c ' of approximately 20° R resuited when the 

thickness ratio increased from 0.02 to 0.10 for the critical point 
behind the shoulder at a free-stream Mach number of 1.00. A decrease in 
T 0 ,c approximately 14.5° R resulted when the thickness ratio 

increased from 0.02 to 0.10 for the critical point ahead of the shoulder 
at a free-stream Mach number of 1.30. 

Experimental verification of the results presented herein has been 
obtained from data taken in a 4- by 10-inch duct tunnel at a free-stream 
Mach number of 1.36 with a 6 • 82 -percent-thick symmetrical-diamond airfoil 
of 5 -50-inch chord at zero angle of attack. The icing limit was deter- 
mined by successively raising the free-stream static temperature until 
no ice formed on the surface. The results of these tests are summarized 
in the following' table and compared with results calculated by the 
analytical method presented herein: 


Free-stream 

static 

Minimum free-stream 
corresponding to ice 

static temperature 
-free condition, Tq >c 

pressure 

Experimental 

Analv 

tical 

*>o 

(in. Hg) 

Front • Back 

surface surface 

Front 

surface 

Back 

surface 

9.13 

404 

403 

— 

9.41 

422 

— 

422 

9.76 

406 

401 

— 

9.91 

418 

— 

422 

3.91 



447 

— 



10 


NACA IN 2861 


The excellent agreement obtained in several cases is probably somewhat 
fortuitous since the experimental error was estimated to be ±10° F. 

These results , however, show that the calculated icing limit is very 
closely approximated experimentally and, hence, provides considerable 
credence for the calculations presented herein. The run-back icing on 
the rear surface, which was indicated by the analysis, was also observed 
experimentally. Visual' observations of the model showed that the surface 
water film flowed smoothly around the shoulder and quickly froze. The 
results presented in the table show that the free-stream static tempera- 
tin^ at which the run -back, freezing occurred experimentally is in good 
agreement with the analytical values. 

The potential hazard to aircraft and missile operation in the tran- 
sonic region may be seen by applying the results previously given for 
the critical points at the shoulder of a diamond airfoil to operation 
under atmospheric conditions likely to be encountered’ If a particular 
relation of pressure and temperature with altitude in the atmosphere is 
assumed, it is possible to calculate the flight Mach number necessary 
for icing protection for a given airfoil thickness as a function of 
altitude. The Mach number necessary for protection for the critical 
point ahead of the shoulder for a 6 -percent-thick diamond airfoil is 
shown in figure 4 for NACA standard atmosphere (reference 11.) . For these 
conditions the region of icing hazard exists from approximately 7500 to 
27,300 feet. If a more seyere temperature relation with altitude is used 
such as might be encountered during winter, that is, a 32° F sea -level 
surface temperature ■ and a lapse rate of 3° F per thousand feet, then the 
curves shown result for the critical points at the shoulder. For these 
conditions the region of hazard exists from sea level to 24,000 feet. 

A Mq of 1.29 is needed at the peak altitude in order to provide 
protection. 

Although the results presented in this report apply only to the 
symmetrical-diamond airfoil at zero angle of attack, they may be used 
as an indication of the region of icing hazard and hence serve as a 
guide for determining where de-icing systems are needed or what flight 
plan may be used to prevent icing. 


' SUMMARY OF RESULTS' 

The results of the analysis of icing limit presented for the 
diamond airfoil at zeroangle of attack in the transonic -supersonic 
speed range are as follows: 

1. The critical point for the whole airfoil in terms of the icing 
limit was located immediately behind the shoulder. The critical point 
on the front surface was located immediately ahead of the shoulder. 

The results indicate that icing may occur at the critical point behind 
the shoulder for Mach numbers as high as 1.4. 
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2. At a constant free -stream Mach number and airfoil thickness, the 
min imum free-stream static temperature corresponding to an ice-free 
condition increases with increasing altitude, 

3. For the. point, immediately ahead of the shoulder the minimum 
free-stream temperature corresponding to an ice -free condition remains 
the same or decreases with increasing thickness ratio at a constant 
value of Mach number and altitude. For the point immediately down- 
stream of the shoulder the free-stream temperature corresponding to 

an ice-free condition increases with increasing thickness at a constant 
Mach number and altitude. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 24, 1952 
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Figure 1. - Variation of free-stream static temperature corresponding to ice-free condition with 

parameter T./T n . (l + 0.176 NL 2 ) for various static -pres sure ratios. Pressure altitude, 

1' 0 , c l 

15,000 feet. 
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Figure 4. - Variation of flight Mach number required for icing protection as 
function of altitude for several assumed atmospheric temperature- pres sure 
relations. Airfoil thickness, 6 percent. 
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ANALYTICAL INVESTIGATION OF ICING LIMIT 
FOR DIAMOND -SHAPED AIRFOIL IN TRANSONIC 
AND SUPERSONIC FLOW. Edmund E. Callaghan 
and John S. Serafini. January 1953. 18p. diagrs. 
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diamond airfoil at zero angle of attack in terms of 
the stream Mach number, stream temperature, and 
pressure altitude. The icing limit is defined as a 
wetted-surface temperature of 32° F an$ is related 
to the stream conditions by the method of Hardy. 

The results show that the point most likely to ice on 
the airfoil lies immediately behind the shoulder and 
is subject to possible icing at Mach numbers as high 
as 1. 4. 
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wetted-surface temperature of 32° F and is related 
to the stream conditions by the method of Hardy. 

The results show that the point most likely to ice on 
the airfoil lies immediately behind the shoulder and 
is subject to possible icing at Mach numbers as high 
as 1. 4. 
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